Enhancement by polydispersity of the biaxial nematic phase 
in a mixture of hard rods and plates 
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The phase diagram of a polydisperse mixture of uniaxial rod-like and plate-like hard paral- 
lelepipeds is determined for aspect ratios k = 5 and 15. All particles have equal volume and 
polydispersity is introduced in a highly symmetric way. The corresponding binary mixture is known 
to have a biaxial phase for k, = 15, but to be unstable against demixing into two uniaxial nematics 
for k — 5. We find that the phase diagram for n = 15 is qualitatively similar to that of the binary 
mixture, regardless the amount of polydispersity, while for K = 5a sufficient amount of polydisper- 
sity stabilizes the biaxial phase. This provides some clues for the design of an experiment in which 
this long searched biaxial phase could be observed. 
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Systems of anisotropic molecules with two symmetry 
axes may form a biaxial nematic phase. In this phase 
molecules align preferentially along two perpendicular 
axes. The biaxial phase is experimentally difficult to 
observe because in systems with biaxial molecules it is 
preempted by smectic or solid phases. This led Alben |l| 
to propose an alternative system which should behave 
similarly, but in which spatial ordering is difficulted: a 
mixture of hard rods and plates. His analysis of this 
system with a mean-field lattice model showed a phase 
diagram with four phases: isotropic fluid (I), rod-like ne- 
matic (N + ), plate-like nematic (N~), and biaxial (B). 
The latter separates the two nematics with second-order 
transition lines when composition is varied from rod-rich 
to plate-rich. Upon increasing concentration, for a rod- 
rich (plate-rich) composition the system first undergoes 
a first-order I-N+ (T-N - ) transition and then a second 
order N + -B (N~-B) transition. At the crossover there is 
simply a continuous I-B transition. Two continuous and 
two first-order transitions meet at this multicritical point. 
At the conclusions of his work, Alben mentions that a 
N + N~ phase separation might replace the B phase, but 
does not considers this possibility in his analysis. A sim- 
ilar phase behavior has later been obtained by other au- 
thors using different models ^, |[ ||, ||. They take into 
account the effect of having free (rather than restricted to 
a lattice) rotations and/or translations. From them one 
learns three main things: first of all, that Alben's phase 
diagram is qualitatively correct; secondly, that there is 
a symmetric mixture, namely that with rods and plates 
having the same molecular volume (hence parallel or per- 
pendicular like particles have the same excluded volume) , 
for which the multicritical point appears more or less at 
cquimolar composition; and thirdly, that the phase dia- 
gram is perfectly symmetric about equimolarity only if 
virial coefficients higher than the second are neglected. 

The first experimental observation of a B phase was 



obtained by Yu and Saupe H in a ternary lyotropic mix- 
ture of potassium laureate, 1-decanol and water. This 
system forms lamellar and cylindrical micelles. By vary- 
ing composition and concentration and B phases ap- 
pear, in a configuration that reproduces Alben's phase 
diagram around the multicritical point. Although this 
was first considered an experimental realization of Al- 
ben's model, it was latter recognized that micelles can 
really change shape from rod-like to plate-like through 
biaxial forms as we move in the phase diagram. A Lan- 
dau theory for a system of shape-changing micelles Q| 
does in fact reproduce even the most peculiar features of 
Yu-Saupe's system (like reentrance in the isotropic phase 
upon increasing concentration) . So the experimental ob- 
servation of the B phase in a mixture of rods and plates 
remains a open problem. 

At the same time van Roij and Mulder || introduced 
an new important element in play. They considered the 
rod-plate mixture version of Zwanzig's model (uniaxial 
parallelepipeds), as well as an expansion of the free en- 
ergy up to the second virial coefficient. Then they reanal- 
ysed the phase diagram with respect to N + -N~ demixing 
for the symmetric mixture, for which the excluded vol- 
ume between unlike particles is minimum (hence phase 
separation is least favored). What they found is that N + - 
N~ phase separation is more stable than the B phase up 
to aspect ratios (long-to-short axis ratio) k w 8.8. Above 
that threshold the phase diagram is like that predicted 
by Alben up to a certain concentration, where phase sep- 
aration again replaces the B phase. The window of sta- 
bility of the B phase is relatively narrow. The driving 
mechanism behind this phase behavior is the larger ex- 
cluded volume between unlike particles in the B phase 
as compared to that between like particles (the rod-plate 
excluded volume divided by the rod-rod one scales as 
k 2 / 3 for large aspect ratios k Q). When the gain in 
free volume compensates the loss in mixing entropy (and 
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this strongly depends on concentration and composition) 
phase separation occurs. This phase behavior was later 
confirmed in simulations of symmetric mixtures of pro- 
late and oblate ellipsoids Q, the only difference being 
the logical asymmetry in the phase diagram of the latter 
(an important difference, though, because it gives rise to 
N~-B coexistence). 

There is a new recent experiment, this time performed 
on a true colloidal rod-plate mixture |^0| . Rod (plate) 
aspect ratio is n kj 10 (15) while plate-to-rod volume 
ratio is 13:1. The mixture is therefore far from being 
symmetric. This system shows isotropic (I) and uniaxial 
nematic phases, as well as biphasic (I-N 1 * 1 ) and triphasic 
(I-N+-N - ) coexistences, but no B phase at all (there ap- 
pear more phases at higher concentrations, but they are 
spatially inhomogeneous |ll|). This picture is consistent 
with the large volume difference between rods and plates, 
and can in fact be explained with a Parsons-Lee density 
functional approximation JllJ. This experiment intro- 
duces, however, another element of unpredicted effect on 
the system: polydispersity. Both rods and plates have 
about 20-30% polydispersity in their axis lengths. For 
this particular system, as theory shows jl]J , this does not 
seem to have any observable qualitative effect other than 
allowing for more than three phase coexistence (which 
is forbidden by Gibbs's phase rule in a two-component 
system) at high concentrations. But as we aim to show 
in this letter, polydispersity has a more drastic effect in 
the phase behavior of the symmetric rod-plate mixture. 

To this purpose, we have extended van Roij and Mul- 
der's model in two ways: (i) including polydispersity 
in a highly symmetric way, in order to minimize triv- 
ial excluded volume effects; and (ii) using fundamental- 
measure theory Jl^, [l3| to model its free energy, which 
for homogeneous phases is equivalent to employing a y- 
expansion exact up to the third virial term (hence the 
phase diagram is expected to be asymmetric). Thus our 
system consists of a mixture of parallelepipeds of size 
L x D x D, all of them with equal volume: LD 2 = 1. If 
we characterize anisotropy by L/D = A, then L = A 2 / 3 
and D = A -1 / 3 . Rods have 1 < A < oo while plates have 
< A < 1. The composition of the mixture is then given 
by a (parent) probability density 



p(A) = A- 1 [C/(A/ K ) + (l-C)/(A K )], 
f(z) = Koia)- 1 exp[-(a/2)(^ 2 + z~% 



(1) 



where K v (a) (a > 0) is the ^th-order modified Bessel 
function, and k > 1. Function f(z) is peaked around 
z sw 1, and it is wider the smaller is a; so /(A/re) is 
peaked around A = k and f(Xn) is peaked around A = 
K . The former represents a polydisperse distribution 
of rods and the latter a distribution of plates, both of 
aspect ratio n. Parameter < £ < 1 allows to tune 
the overall composition of the mixture, since the molar 
fraction of the rods is given by x r (Q) = p(A) dX (and 



that of plates by x p (() = 1 — x r (Q, of course). The 
probability density p(X) is chosen so that if dg(X, C) = 
p(X)dX, then \dg(X,()\ = \dgiX~ 1 ,1 — ()\; in other words, 
the fraction of rods and plates within the same interval 
of aspect ratios are in the proportion x r (() : x p (l — £). 
The moments of this distribution are given by (A m ) = 
K m/2 (a)iY (a)- 1 [C« m + (1 - 0>^ m ], explicitly showing 
the symmetry of the mixture. 

Fundamental-measure approximation for a multicom- 
ponent mixture of hard uniaxial parallelepipeds amounts 
to taking for the excess (over the ideal) free energy den- 
sity (in kT units) g|, 
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pf denoting the number density of parallelepipeds of 
species i with their symmetry axes oriented along the 
direction /i(= x, y, z). Specializing for our system, where 
species are labeled by the continuous parameter A, 
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where p(A) = P^W- The total free energy density is 
then given by 
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Details for determining phase equilibria in a system 
like this have been reported elsewhere |lj, [Tq] . In brief, 
if at a given p there is only one phase present, its equilib- 
rium composition is determined by minimizing §) with 
respect to p M (A) under the constraint p(A) = pp{X), with 
p(A) given by (|l]). This amounts to solving the system of 
equations (p = x, y, z) 



p"(A) = pp(X) 
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If to > 1 phases are present in mutual equilibrium, then 
labeling them with a(= 1, . . . , to) their density distribu- 
tions must verify the 'lever rule' 



(5) 



v a denoting the fraction of volume ocupied by the ath 
phase (Ea v a = !)■ Chemical equilibrium then yields 



p£(A) = pp(A) 



Eb^bE^e 



(6) 
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and this system of equations is completed by the equality 
of osmotic pressures between every pair of phases, where 
osmotic pressure is obtained from (0) as 



pu = 



1 - 9 ' (1 " 9? 
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First of all, as this theory has contributions of virial 
terms higher than the second, we have checked van Roij 
and Mulder's results |p[. So we have first chosen a bidis- 
perse mixture (pure rods and pure plates). The resulting 
phase diagrams are qualitatively the same, apart from 
some asymmetry induced by those higher virials: for 
k = 15 there is a small B window bounded above by N + - 
N~ demixing, whereas no B phase appears for k = 5. 
In both cases, the multicritical point is found very near 
x r = 1/2. An important difference produced by the 
asymmetry is the replacement of the very narrow region 
of triphasic coexistence found by van Roij and Mulder 
by a simple N + -B coexistence. It is also worth mention- 
ing that the asymmetry we obtain mirrors that of the 
simulations of prolate and oblate ellipsoids (|] (they ob- 
serve N~-B coexistence instead). This means that the 
asymmetry must be strongly influenced by the details of 
the model (mainly the shape of the molecules and the 
restriction of orientations). 

We have next made the mixture polydisperse. The 
breadth of the two peaks in p(X) is controlled by a in (Q). 
A quantitative characterization of the polydispersity can 
be given if we determine the dispersion in L and D as 
obtained from p(A) for £ = or 1. This yields 




K v /s{a)K Q {a) _ 



(8) 



where v = 2 for and vd = 1 for Ad- For k = 15 we 
have chosen a = 1 (A L = 0.288, Ad = 0.143). The re- 
sulting phase diagram (p vs. x r ) is plotted in Fig. |l|. The 
global picture (Fig. |l]a) resembles very much the bidis- 
perse case: there are first order I-N ± transitions; the two 
uniaxial nematics are separated by a B phase through 
second order transition lines; roughly above a threshold 
density the latter is replaced by a wide region of N + -N~ 
coexistence; and there is a multicritical point where the 
two I-N* first order transitions and the two N ± -B sec- 
ond order ones meet (very close to x r = 1/2). There is 
a difference with respect to the bidisperse phase diagram 
in that the lines delimiting coexistence regions do not 
correspond to coexisting states, but are the cloud lines 
(defined by points where a incipient new phase is form- 
ing) of the corresponding coexistences. For the sake of 
clarity, the shadow lines (density and composition of the 
incipient phase) are not represented. Notice that p is the 
global density of the system, and does not coincide with 
that any of the coexisting phases (except at the cloud 
lines of biphasic coexistence). The details (Figs. |Tjt> and 
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FIG. 1: Phase diagram (global density, p, vs. rod molar frac- 
tion, x r ) and two details of a polydisperse rod-plate mix- 
ture with length and breadth polydispersity Al = 0.288 and 
Ad = 0.143, and aspect ratio k = 15. Phases are labelled I 
(isotropic), N + (rod- like nematic), N~ (plate-like nematic) 
and B (biaxial nematic). Dotted lines mark second-order 
phase transitions and full lines delimite coexistence regions. 



c) show that right above the B phase there are regions of 
B-N ± as well as triphasic coexistence. Also remarkable 
are the two second-order transitions separating the B-N* 
coexistence regions from the N + -N _ one, where upon in- 
creasing concentration the biaxial order parameter van- 
ishes continuously, hence transforming the B phase into 
the second uniaxial nematic. These lines have no ana- 
logue in the bidisperse system. Finally, a global effect of 
polydispersity is to lower the phase diagram in densities 
and to make it more symmetric (probably two related ef- 
fects, because the lower the density the less relevant the 
higher virial terms). 

We have checked that the region of stability of the B 
phase is not appreciably affected by an increase in poly- 
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FIG. 2: Same as Fig. [l] but for aspect ratio k = 5 and length 
and breath polydispersity Al = 0.610 and Ad = 0.302. 



dispersity (using a = 0.1, i.e. A L = 0.610, A D = 0.302), 
so polydispersity does not seem to inhibit B ordering. In 
fact, it acts otherwise, favouring B ordering against N + - 
N~ demixing, as the case k = 5 illustrates (see Fig. ||). 
Again the phase diagram of this case looks very much 
like that of the bidisperse case if we take a = 1, i.e. 
no B phase appears. However, increasing polydispersity 
up to a — 0.1 we again obtain a region where B order- 
ing is more stable than N + -N~ demixing (Fig. ||a). The 
B phase is limited from above by a region of B-N~ co- 
existence, which upon increasing density again becomes 
N + -N _ through a second order transition. 

This enhancement of biaxial ordering is the most srik- 
ing effect of polydispersity. To understand why it is so 
we must find a mechanism by which B ordering is en- 
tropically favoured w.r.t. N + -N _ demixing. This mech- 
anism is two-folded: on the one hand mixing entropy 
increases upon increasing polydispersity, thus penalizing 
demixing; but on the other hand, polydispersity increases 
the gain in free volume of the B phase w.r.t. that of 
uniaxial nematics (the average rod-plate excluded vol- 
ume divided by the average rod-rod one in a perfect B 
phase scales as K 2 / 3 e _cA , c being a positive constant, 
for large aspect ratios k and polydispersities in the range 
1 < &(L,D) < VIok)- 

This effect can be exploited in the design of an exper- 
iment to observe the B phase. The fabrication of rod- 
like and plate-like colloidal particles could follow simi- 
lar procedures to those employed by van der Kooij and 
Lekkerkerker (T^]. Then, the mean size and polydisper- 
sity of the rods and the plates can be controlled by in- 



ducing successive fractionations (e.g. by adding a non- 
adsorbing polymer ]l6|), starting off from the appropri- 
ate parent distributions so as to obtain the desired final 
values. The key is to produce as symmetric (in parti- 
cle volume) polydisperse mixtures as possible, since this 
seems a crucial ingredient for the stability of the B phase, 
with or without polydispersity. How destabilizing are the 
asymmetry of the mixture or the existene of a particle- 
volume distribution (motivated, for instance, by having 
independent polydispersities in length and thickness) is 
something that has yet to be quantified, but we expect 
these results to be robust against perturbations in these 
two directions. 

A final point to be addressed is that we have not taken 
into account possible nonuniform phases in the determi- 
nation of the phase diagram. According to the availble 
experiments |Io| l these phases appear at sufficiently high 
concentration, while the relevant parts of the phase dia- 
grams we present all occur at rather low concentrations. 
So it seems unlikely that these parts are affected by the 
presence of inhomogeneous phases, but this has yet to be 
analyzed with some care. 
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